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Abstract 

By examining BHDs from the general quantum-field-theoretical point of view we derive a num- 
ber of generalizations of the standard analysis of their response. In particular, by allowing for 
interactions restricted in time (by a smooth function) we present general expressions for BHD's 
output (in which the usual, simplifying limits can but need not to be taken). Moreover we point 
out the need for non strictly-monochromatic local oscillators (i.e. the need for "pulsed" ones) in 
order to have well-defined QFT-observables (the products of which, eg. have finite vacuum ex- 
pectation values). Furthermore we show how the analysis of the detectors generalize to situations 
with BHDs in waveguides, Casimir cavities, or other time-independent but inhomogeneous e and 
li. The general treatment also allows us to comment on important QFT features of the detector- 
observables such as locality (i.e. commuting for causally separated measurements) and non-null 
vacuum response. This leads to the conclusion that balanced homodyne-type detectors (and not 
single-photodiode-like ones) are the appropriate tools in testing intriguing QFT predictions like 
negative (sub-vacuum) energy-densities. Finally by recalling results on field-autocorrelation func- 
tions for QFT in Casimir cavities we show that interesting effects (large reductions of fluctuations) 
are to be expected if a version of BHDs were to be placed in such cavities. 
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I. INTRODUCTION 



In quantum field theory the vacuum is not as empty as in the classical physics. The 
vacuum expectation values of positive operators, for instance of the square of the smeared 
electric field, are non-vanishing and there are plenty of states with the expectation values 
lower then the vacuum ones. These states are locally "darker than vacuum" . While the un- 
derstanding of this phenomenon progressed constantly over the last decades, with a number 
of surprising results having been discovered [20], attempts were also made to see this non- 
classical behavior in experiments. Naturally, in order to see sub- vacuum fluctuations it is at 
least necessary to have detectors capable of seeing and quantifying the vacuum fluctuations. 
Such detectors, known as balanced homodyne detectors (BHD) with local oscillators, were 
first proposed in the context of quantum optics by Chen and Yuan [2] . With their help the 
so-called squeezed states of light are seen to exhibit regions with sub-vacuum electric field 
fluctuations (eg. [3]). 

In this paper we reconstruct the quantum-field-theoretical observables corresponding 
to photodiodes and balanced homodyne detectors with local oscillators and discuss their 
quantum-field-theoretical properties. The paper is organized as follows: in the part A of the 
second section we investigate the interaction of a simple quantum system with the quantum 
electromagnetic field and derive a formula for the response of a photodiode. The part B 
deals with an arrangement of two photodiodes known in quantum optics as the balanced 
homodyne detector. This part contains four subsections dealing with generalizations of the 
standard derivation of the observable corresponding to measurements on BHDs. In partic- 
ular we show in which asymptotic sense the expectation values of the output of BHDs, and 
its variance, are related to the one- and two-point functions of the quantum electromagnetic 
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field (subsections 1 and 4). We generalize these results to the case of non- monochromatic 
local oscillators (subsection 2) and to the case of BHDs placed in inhomogeneous environ- 
ments (waveguides, Casimir cavities, inhomogeneous media; section 4). The utility of these 
generalizations is shown in the appendix A, where we analyze the expected response of BHDs 
placed in Casimir cavities (consisting of two parallel, perfectly conducting plates), and argue 
that substantial, time- independent, reductions of fluctuations (aka sub- vacuum expectation 
values of the square of the electric field operator) are expected. The third section contains a 
number of comments on the nature of BHD-measurements and provides a relation to more 
abstract QFT results (vacuum response, causality and the Reeh-Schlieder theorem). 

II. DETECTORS 
A. Photodiode 

Let us start the discussion with a photodiode, on which a photodetection process takes 
place in a semiconductor region (usually the I region of a PIN diode). If light is shed upon 
this region, electron-hole pairs are created. Upon creation, these charges propagate towards 
and eventually reach the opposite sides of the diode, where they are collected. In this paper 
we shall model the detection as an excitation of an electron, initially in a relatively well- 
localized bound state ip (eg. of an isolated donor atom), to the continuum of excited states, 

We will adapt here the of standard treatment of the detection process, presented eg. in 
the complement A n of [4], or in the section 4.2 of [5]. Let the dynamics of the electron be 
generated by the Hamiltonian 

H = H + V mt (1) 
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with 



H = p 2 /2m* + 1/(x) 



(2) 



ipo an d V'(q) being eigenstates of if > an d the dip ole- approximation interaction-Hamiltonian 



where ex 1 = d l is the dipole-moment operator and Ei(t, x) is the (quantum) electric-field 
operator. We shall assume, that this field-operator can be replaced by the electric-field 
operator taken at the point, x, where the donor-atom resides, 



This simplifying assumption will be argued unnecessary for the BHD in section III B. The 
evolution will be determined from the first order time-dependent perturbation theory. We 
shall use the formulation, [6, 7], in which the interaction Hamiltonian is multiplied with a 
smooth function of time, g(t), equal to one in the time interval during which the measurement 
takes place[21], and vanishes rapidly just outside of this interval. The probability of exciting 
the electron by a given initial state S of the quantum field is just the expectation value of 
the (time-evolved) projection operator on the excited states of the electron, say P e ® 1, taken 
with respect to the state tp Q ® S of the full system. The well-known first order perturbative 



V int = ex l Ei(t,^), 



(3) 




(4) 
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expression for this probability is 



PD(<?,x) = (^o ® s\u; (P e 1) U g \4> ® S) = 

/oo /■ 
^r^tfedr / rf/i(q) ty„, e^ s tf e^>(q)) ty(q), <P e^^o) (^(r,x)^(s, 
-oo «/ 

where C/ g is the evolution operator with the interaction controlled by g(t), and the measure 
d/i(q), specifies the final, scattering states of the electron. The above expression can be 
written in a compact form: 

PD((?,x) = J A^(q)e^ (r - s) ^(qT^(q) (^(r, x)Ej(s, x)> 5 , (6) 

where symbolically d 2 g = g(r)g(s) dsdr, while w q = -E'(q) — -Eo stands for the unperturbed 
(w.r.t. i? ) energy difference between the excited and the ground state of the electron and 
<# (q) = (^(q), exi i])q). It is convenient to introduce a symbolic star notation, such that 

/oo 
d 2 g (£(r,x)*£( S ,x)) s (7) 
-oo 

B. Balanced homodyne detector (BHD) with a local oscillator 

The primary purpose of a BHD with a local oscillator (LO) is to investigate the states 
of the radiation field, which alone hardly trigger photodiode's response, eg. vacuum/ground 
states. Consider the setup presented on the figure 1. The observable corresponding to the 
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charge collected at V is the difference of the outputs of the photodiodes: 



J = PD(<7,x)-PD(<7,y). (8) 




PBS1 HWP PBS2 



FIG. 1: Balanced homodyne detector with a local oscillator. The linearly polarized signal field, S, 
(if present) is "blended" with a coherent state (LO), which is polarized orthogonally to S, on the 
polarizing beam splitter (PBS1). The half wave plate (HWP), reflects the planes of polarization 
with respect to its optical axis, thereby inducing 7r/4 shift of the plane of polarization of the signal 
field S. The subsequent PBS2 separates the two orthogonally polarized signals, which are detected 
at the photodiodes PDx and PDy. The charge collected at V provides a measure of the expectation 
value of the observable J (and of its higher moments). Note, that the setup is arranged in such a 
way, that if S happened to be a monochromatic coherent state, it would be phase-matched to the 
LO at the point x, but shifted in phase by it at the point y. 

In experiments with BHD the state under investigation is "blended" with a strong, pre- 
cisely controlled coherent state (LO). Let S denote the state of the radiation field under 
investigation (that is: the state of the quantum field without the local oscillator present). In 
a typical experimental setup (fig. 1) the polarising beam splitter, PBS1, superposes phys- 
ically the state S and the LO. The blending has a "coherent character" and therefore we 
adapt one of the definitions of the coherent state to describe it: 



(P[Ei(t,x)]) (StF) = (P[Fi(t,x) + E i (t,x)]) s , 
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where 



*i(t,x) = f (E t (t,x)) F (9) 



denotes the electric field of the LO. The symbol P above stands for an arbitrary observable 
(polynomial) constructed out of the (smeared) electric field operators, while (S, F) denotes 
the state resulting from the blending of the local oscillator coherent field into the state S 
under consideration. We choose the LO field to be monochromatic and linearly polarized, 
so that at every point x we may write 

Fi(r, x) = Ki [e-^-' ) + e Mr-to)] ? (10) 

with a real vector K { and the frequency u > 0. The phase t , which depends on the point 
x, can be varied at will in experiments. Per definition the detector is arranged (balanced) 
is such a way that [22] 

F l (r,x) = -F,(r,y). (11) 
The expectation value of J for the detector balanced in this way is 

(J) {s ,f) = (J)s + J d 2 g {[F(r,x) * (£(s,x))s + <£(r,x)>s* F(s,x)]- 

- [F(T,y)*{E(s,y)) s + {E(T,y)) s *F(s,y)]} (12) 

(the terms proportional to F 2 cancel.) The quantum electric field is 

^(s,x) = ~d s Ms, x) = -d s / 4= [e- tps a t (p) + e^a*(p)} , (13) 
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where p — |p| and 

Mp), a*(k)] = 8(p - k) - (14) 
It can also be written briefly (by defining appropriate &«(p)): 

^(s,x) = J d 3 p [e- ip %(p) + e ips b*(p)} . (15) 

Let us now symbolically perform the r and s integrations in the operators of the form 
J d 2 gF -k E, the expectation values of which appear in eq. (12). The integration leads to 
terms involving 

(?(^ q =p u) g (TP - w q ) • (16) 

Of these, clearly the term g(u; q — oj)g{p — cj q ) dominates the other three, because u,Uq,p 
are all greater than zero, and the Fourier transform of g(t) is concentrated around the 
zero frequency with a rapid decay away from it. This, as explained in the appendix A, 
essentially enforces the restrictions c<j q = uj and p = uj in the integrations over d/x(q) and 
d 3 p respectively. 

1. Asymptotic response of BHDs 

At this point let us make a simplifying assumption that we are dealing with an electron 
in a spherically symmetric initial state, ip = ^ (l x l)i with the energy E which after ex- 
citation becomes a plane wave ^ q = exp(iqx), and thus cfyt(q) = d 3 q and a fixed value of 
c<j q = q 2 /2m* — E translates into a fixed value of q = |q|. In a real semiconductor both 
of the functions i/jq, ip(q) are different but, of course, for a given crystal they are very well 
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known. With these assumptions we find 



d*(q) = eq%q), (17) 

with a certain function h(q). The d 3 q integration is now straightforward; with the integral 
over the directions, / dfl q l q j = ^-S ij q 2 , we finally get for the leading contribution from 
f(PgF*E: 

J d 2 gF(r)*E(s) = __ m ^|/,( go )| 2 ^ J d 3 P 5(p-u J )b*(p)e^ (18) 
where q = ^/2m*(w + E ). This expression is equivalent to 



J d 2 gF(T)*E(s)=A(iu) K'E'ito) 



with E^~(to)\ u meaning the negative frequency part of the electric field operator restricted to 
the frequency iv, taken at the time to, and with the abbreviation A{u) = ^^m* ql\h{qo)\ 2 . 
Analogously we find 

j d 2 g E(r) * F(s) = A(u) K l E+(t )\ u . 

The expectation value of the output of the BHD detector, (J)(s,f), eq. (12), is therefore 
given to the leading order by 

{J) (s ,F) = A(u)K i <^(t ,x)U + ^(to,y)U> 5 . (19) 
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2. Limit T — > oo and monochromatic LO 

It must be stressed, that the limit T — > oo taken here with the assumption of monochro- 
maticity of the local oscillator field, which leads to (18), is not physical and mathematically 
troublesome (intergating the response for infinitely long times / smearing with singular test 
functions leading to products of delta functions for products of operators of the type (18)). 
To avoid these troubles one should replace K % by a frequency- dependent function ki(u) and 
integrate (18) over uo. The function k^u) can be for instance such as to correspond to 
"pulsed" almost-monochromatic local oscillators; the equation (18) is then replaced by 



with a smooth k l {p) sharply peaked around the frequency u of the local oscillator. 

3. Generalization to non-trivial geometries or inhomogeneous media 

For quantum fields in non-trivial static environments (waveguides, cavities or materials 
with position-dependent e(x) and /i(x)), the electric-field operators possess a generalization 
of the decomposition (15) 



where the multi-index p a corresponds to the parameters of electromagnetic waves supported 
by the environment (such as types of waves and their wave vectors), ?/>j(p a ,x) denotes the 
(electric field of the) solution of the Maxwell equations with these parameters and du(p a ) 
denotes the appropriate measure. The frequency, u>(p a ), is given by an appropriate dispersion 




(20) 




dis(p a ) e- iu{pa)t b{p a )^{p a + e iu}{pa)t b*(p a )A(p a ,x) , 



(21) 
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relation. The functions i/)i(/,x) are normalized so that [b(p a ),b*(k a )] = 5(p a — k a ). The 
generalization of eq. (20) is then 

J d 2 gF(r)*E(s) = A{uo) J du(p a ) k l \oo (p a )]b* {p a )^{p a , x) = A(u)E-(t )\ k(u) , (22) 

which is the negative-frequency part of the field operator restricted in frequencies and po- 
larizations by k l {uj). This seemingly straightforward generalization has the profound effect 
of making the commutator of positive and negative frequency part of the field operators po- 
sition (x) dependent. As a consequence there may be reduced fluctuations of both so-called 
quadratures of the quantum field. An example of this is presented in the appendix A. 

4- Variance of the BHDs output 

For many interesting states S, the expectation value of the electric field operator (and 
thus the leading term of (J)(s,f)) vanishes, for instance if S is the ground/vacuum state, or 
an eigenstate of the photon-number operator, or a superposition of pairs of photons - eg. the 
squeezed state. Let us present the expectation value of the square of J in the ideal case of 
monochromatic LO and homogeneous media. The calculation is trivial as we manipulate the 
operators, under the expectation value, which are a product of two independent expressions. 
The leading term of (J 2 ) is given by 

(J 2 ) { s,F) = A\u)K i K^[E t (to,x)l + E t (to,y)l][E^ (23) 
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III. DISCUSSION 



A. Properties of (J 2 ) 
We note the following: 

• (J 2 ) is proportional to K 2 , that is to the power of the local oscillator field; by taking 
measurements with growing values of K 2 the value of the fluctuations of the electric 
field in the state S can be quantitatively estimated. This is the standard signature 
checked for in experiments, see eg. [3, 9]. 

• Should one wish to compare (J 2 ) of two different states on the decibel scale, the 
dependence on the yet unspecified parameter, A, and LO power drops out, and the 
result depends only on the (frequency-restricted) two-point functions of the states 
under consideration. 

• If the state S can be assumed to have certain symmetries, a further simplification of 
the form of (J 2 ) is possible. For instance if S = fl is the Poincare invariant (usual) 
vacuum, then 

(J 2 ) is , F) = 2A 2 H^^'(^(t ,x)U^(t ,x)U) s (24) 

(The spatial orientation of the photodiodes plays a role here.) 

• All the operators are restricted to the frequency, u, of the LO; therefore time is here 
2ty/uj periodic, and thus t G [0, 2ti/uj) is the measure of time with 2tt/uj corresponding 
to the period of the LO wave. 
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B. Interaction Hamiltonians 



With regard to the interaction Hamiltonian which has been employed, (3), we note that 
we could have taken it without the restriction on the electric- field operator, eq. (4); its 
frequencies would nonetheless get restricted to the frequency of LO (at the BHD) and the 
large wavelength of light waves of this frequency (against the essential support of ip ) would 
lead to the same formulas for the output of the BHD. Moreover, because of the standard 
relation 

(V>(q), P 3 ' V'o) = Mn*w q • (V>(q), x 3 ip ), 

which is generally valid if tp and ip(q) are eigenstates of the Hamiltonian H , the same 
expressions for the expectation values of J and J 2 would have been obtained had we started 
with the interaction Hamiltonian of the form -^A l (t, x) pi. 

C. Note on the locality of measurements 

Let us assume that ?/>o( x ) is supported in a compact region [23] and that g(t) is also of 
a compact support (in time). Then the quantum-field-theoretical observable correspond- 
ing to PD(g,x) is just an integral (over dfi(q)) of a product of two electric field operators 
smeared with test functions supported in the spacetime region supp i[) x supp g. Such an 
observable is QFT-localized in the smallest double-cone [24] containing this region. As a con- 
sequence, two such observables localized in causally separated double-cones will commute 
(and thus have their probabilistic distributions unconstrained by Heisenberg-like uncertainty 
relations). Moreover, since the operator corresponding to PD(g,x) is clearly positive, it can- 
not annihilate the vacuum because otherwise by the theorem of Reeh and Schlieder[10, 11] it 
would necessarily be a null operator (vacuum is separating for the local observables). Thus, 
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strictly, photodiodes are not "photon counters" , as they exhibit non-null vacuum response. 
These "vacuum-effects" were also encountered in [8, 12], and are related to causality prob- 
lems of Hagerfeldt [13] (see also the subsequent clarification by Buchholz and Yngvason 
[14])- 

IV. CONCLUSIONS AND OUTLOOK 

We have corroborated the result, that the output of a balanced homodyne detector pro- 
vides a quantitative measure of the n-point functions of the states of the quantum electric 
field. Moreover we have shown how this result can be generalized to the case of detectors 
placed in inhomogeneous environments. Non-trivial effect are expected in such a case, as is 
shown in the appendix A dealing with quantum fields and BHDs outputs in Casimir cavities 
(see also [16-19]). Our investigation supports the view, that BHD-like detectors are suitable 
for measurements of intriguing QFT effect such as negative energy densities present for in- 
stance for squeezed states, and thus provide a tool for testing recent theoretical predictions 
associated with Quantum Energy Inequalities [1, 15]. 
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APPENDIX A: BALANCED HOMODYNE DETECTORS IN CASIMIR SETUPS 



In this appendix we will derive a prediction for the output of a BHD placed between 
Casimir plates. In order to do so we will recall the results of [18] and provide their relation 
to the characteristics of the output of the detector. In ref. [18] autocorrelation-functions 
(two-point-functions) were derived for the ground and thermal state of the quantum electro- 
magnetic field between perfectly conducting parallel plates. Let us take the ground state G 
for simplicity. The autocorrelation function for the electric field operators, (Ei(x)Ej(x')) G 
is obtained (Eq. (2.25c) of [18]) by partially differentiating (w.r.t. the first argument) the 
image-sum-functions 



where L is twice the distance between the plates which are perpendicular to the re-direction. 
Specifically for electric field operators in the y-direction (tangential to the plates) we have 



where differentiations are performed with respect to the first point (x) only. Let us now 
set the points so, that x lie at the same spatial position as x' but by s later in time: 
x = (s, x, 0, 0), x' — (0, x, 0, 0), where x is the distance from the first pate. We find 



oo 



(x t x> - nL) 2 + (y- y') 2 + (z - z') 2 - (t - t') 2 



1 



n=— oo 



(E s (x)E y (x')) G = (dl + S?) [F~(x, x') - F+(x, x' 



)]■ 



1 



oc 



(nL) 2 + s 2 (nL - 2x) 2 + s 2 



(E y (x)E y (x')) G 



TV 



- y 

2 Z—^ 



(s - nL) 3 (s + nL) 3 (s - nL + 2x) 3 (s + nL - 2x) 3 



n=— oo 
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which is a function of s and x. The Fourier transform of this function with respect to s 
gives the frequency spectrum (eq. (2.30) of ref. [18]), 

i poo 3 00 

<tg(w,x) = — / ^e ms (E y (x)E y (x')) G = 4^ 5^ {Q[omL] - Q[c;(raL - 2x)]}, 

— 00 n=— oo 

where the function Q is given by 



^. . SIM cosx SIM 
Q(x) = + 



X 



X' 



X 



3 ' 



We note, that the n = term of the first summand, equal to uj 3 /Qit 2 corresponds to the 
Minkowski- vacuum, fl, spectral density (i.e. to the situation without plates). Hacyan et 
al. have summed up the series and expressed the n-sums in terms of piecewise elementary 
functions (polynomials and trigonometric functions of a;). For the purpose of this paper let 
us merely plot vacuum and Casimir spectral densities for L = 2/im (plate separation lfim). 



-[— ] 

c 1 ' 



FIG. 2: The normalized difference (<tg(w,x) — (Jji(w,x))/(7sj(w,x) between vacuum and Casimir 
spectral densities for u/c G [0, 47r]//m~ 1 . Black color corresponds to negative values (sub- vacuum 
spectral densities). For uoL < 2ir Casimir spectral density vanishes, <tg(w,x) = 0, for all x- 
Discontinuities appear at luL = 2nir. 



Let us finally make a connection to the observables measured by BHD's. For quantum 
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FIG. 3: Frequency-dependence of the suppression of fluctuations in the dB scale, 
10 • log 10 [o"g-(w,x)/cq(u;, x)], for the ground state in the Casimir geometry. Plotted curves corre- 
spond to x = 0.25 fxm (solid) and x = 0.5 fim (dashed). 

fields in non-trivial static environments we find, that (in the notation of section II. B. 3) 

<tg(w,x) = f dv(p a ) du(k a ) 5 [u(p a ) - u] \^ y (p a ,x)\ 2 5(p a - k a ). (Al) 

The expectation value of the output of the BHD in the ground state vanishes, while its 
variance, (23), contains two parts, one of which depends on the relative position of the 
diodes (this part comes from the term (Ei(t,XL)Ej(t,y))). The position- independet part is 
simply related to the spectral density: 



{J )indep = ^ (^Lo)(^(t0,x)UH^(t0,x)|^) +^( t 0,y)|fcH^(t0,y)UH> = ( A2 ) 



= 2A\u LO ) J duja(u,x)\k y (uj)\ 2 , 



(A3) 



where ujlo is the local-oscillator frequency. Thus in the usual situation where k y {ui) is sharply 
peaked around ui^o the spectral densities presented above are proportional to the position- 
independent (and LO-phase, t , independent) output of a balanced homodyne detector. 
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APPENDIX B: PROPERTIES OF THE SMEARING FUNCTIONS 



In the paper, the arbitrary functions of time, g(t), multiplied the interaction Hamiltonian. 
These functions should be set to one in the spacetime region, where the actual interaction 
takes place. Thus we view them as being smooth and essentially equal to one in the interval 
[—T/2, T/2], where T is taken sufficiently large. Let us take for instance 



g (t) = - {tanh[a(t + T/2)] - tanh[a(t - T/2)]} . 



we find easily (using the convention g(uj) = f e lwt g{t) dt), 



tanh(c<;) 



sinh(7rc<j/2) ' 
and 

7r sin(u;T/2) 



a sinh(7ru;/2a) 

This function evidently decays rapidly for u large with respect to a. Moreover g(0) = T, 
and due to the known representation of the Dirac's delta 



sin(^T) 
hm = no (uj) 

T^oo UJ 



we have 

lim g{uj) = 2tt8{ui). 

This property is generic; any positive real function of time, gr(t), which converges weakly 
(on the space of functions of rapid decay) to one as the parameter T — > oo has the property 
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gr{oj) — > 2ix5(uj). If gr(t) is smooth, c/t(uj) will be of rapid decay (for finite T). 
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